UNCLASSIFIED 


AD  NUMBER 


AD809147 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  FEB  1967.  Other  requests  shall 
be  referred  to  Rome  Air  Development 
Center,  Griff iss  AFB,  NY. 


AUTHORITY 


RADC  ltr  17  Sep  1971 


THIS  PAGE  IS  UNCLASSIFIED 


REPUCTO  R  ANTENNA  ZOQ&  TECH  NICHES 

•  .  .  _  •.  Lerpe  K-  OeSfxe 

Peter  A,  Metartes 
,J\  Kj  ■  v  -  Geergfr  E.  Sktffetl) 

1  A?rb«t»e*  Insfrumetits  t-ot^uftry 

v>*  .  ' 


*  <•«  t 

T^CUNJCAt  REPORT  RG,  F^DC-TR-  66-714 

r  •  *  , .  V36?:,. 


ia  Imijeet  to  kpecttfl  '•’ 
.  **!*>«  C  '  Str&ls  *fitj  e«cfe  liasamittil 

"  *0  fee^t»r»br»*nw>6flUi^^^ 
tiesil*  qz  npmtmWv**  tfwafci  tuy 


REFLECTOR  ANTENNA  ZOOM  TECHNIQUES  l 

Lome  K.  DeSize  ; 

Peter  A.Mclnnes  , 

George  E.  Skafiill 

Airborne  ! ns frtmcnfs  Laboratory 


This  document  is  subject  to  special 
export  controls  and  each  transmittal 
to  foreign  governments,  foreign  na¬ 
tionals  or  representatives  thereto  may 
be  made  only  with  prior  approval  of 
RADCtEHLI),  GAFB,  N.Y-  13440 


U* 


1* 


•* 


xnc,  ears,  j.r„  a  x*r  Sr-na. 


i*-.  '  ' 

ff«o  •-  j 


FOREWORD 


This  report  was  prepared  by  the  Airborne  instruments 
Laboratory,  A  Division  of  Cutler-Hammer,  Inc, -  Deer  Park.  N^w  York, 
under  Contract  AF  30(6D2)-4179,  Project  04,  Task  4506,  and'  describes 
the  work  performed  during  the  period  from  22  June  to  22  September  1966. 
The  contractor’s  internal  report  number  is  1031-1-2.  Mr.  Donald  A. 
Hildebrand  is  die  RADC  Project  Engineer  on  this  contract.  The  authors 
report  was  submitted  by  the  authors  October  1966. 

Document  contains  information  erroargoed  from  release  to  Sino- 
Soviet  Bloc  Countries  by  AFR  400-10,  "Strategic  Trade  Control  Program." 

This  report  las  been  reviewed  and  is  approved. 


Approved: 

DONALD  A.  HILDEBRAND 
Project  Engineer 

Antenna  &  Coherent  Optical  Section 


ABSTRACT 


This  report  describes  the  work  performed  during  the 
second  quarter  of  a  study  of  a  reflector  antenna  that  provides  j scorn 
(variable  beamwidih)  and  scan  capability  using  controlled  aperture 
amplitude  and  phase.  The  antenna  consists  of  a  primary  reflector 
(paraboloid)  and  a  secondary  reflector/lens.  It  operates  as  a  lens  in 
conjunction  with  one  feed  for  scanning  in  the  receive  mode  and  as  a 
reflector  in  conjunction  with  another  feed  for  zooming  in  the  transmit 
mode.  Switching  between  a  zooming  transmit  mode  mid  a  scanning 
receive  mode  results  in  a  versatile  radar  antenna  with  an  inherent 
duplexing  capability.  The  performance  of  this  antenna  system  is  be¬ 
ing  analyzed  numerically  with  the  aid  of  digital  computers.  This 
report  describes  the  development  of  a  computer  program  for  solving 
the  Fraunhofer  aperture  integral  efficiently  and  with  good  accuracy. 
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SECTION  I 
INTRODUCTION 


For  several  years  the  Rome  Air  Development  Center  (RADC) 
has  sponsored  *wo  techniques  that  make  use  of  some  form  of  a  reflector 
antenna  in  conjunction  with  aperture  phase  and  amplitude  control  to 
provide  zooming  (variable  beamwidth)  and/or  scanning. 

One  technique,  developed  by  Airborne  Instruments  Labora¬ 
tory  (AIL),  uses  a  cluster  of  feeds  placed  on  a  spherical  surface  con¬ 
centric  with  the  focal  point  of  a  paraboloidal  reflector.  Each  feed  is 
then  independently  controlled  in  phase  and  amplitude  to  provide  zooming 
and/or  scanning  (reference  1).  The  other  technique,  developed  by 
Blass  Antenna  Electronics,  uses  a  flat  reflector  consisting  of  individual 
waveguide  elements  whose  phase  can  be  controlled  by  diode  switches 
that  change  the  position  of  the  short  circuit  in  the  waveguide.  Proper 
programming  of  the  switches  allows  the  beam  from  the  reflector  to  be 
varied  in  width  and/or  scanned  (reference  2).  This  system  is  commonly 
called  a  reflectarray. 

The  purpose  of  this  study  is  to  consider  another  technique, 
different  from  the  two  previously  mentioned,  that  will  use  a  reflector 
and  result  in  a  variable  beamwidth  and  a  scanning  capability.  The 
specific  technique  to  be  investigated  attempts  to  use  the  advantages  of 
the  two  previously  mentioned  systems  in  a  single  system  called  (for 
lack  of  a  better  name)  a  hybrid  system.  This  hybrid  system  uses  the 
optically  fed  phased  array  approach  of  the  Blass  system  together  with 
the  beam  magnification  and  phase  correction  advantages  of  the  AIL  sys¬ 
tem  (Figure  1). 

The  hybrid  system  consists  of  a  primary  reflector  (parabo¬ 
loid)  a  secondary  reflector  acting  as  both  a  phase- controlled  reflector  and 
a  phase-controlled  lens  (shape  to  be  determined),  and  two  feed  systems. 
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VERTEX  FOCAL  POINT  FEED 

FEED 


FIGURE  1.  HYBRID  ANTENNA  SYSTEM 

The  focal-point  feed  system  is  used  in  the  receive  mode  with  the  secon¬ 
dary  z’eflector  acting  as  a  lens;  the  vertex  feed  system  is  used  in  the 
transmit  mode  with  the  secondary  reflector  acting  as  a  reflector.  Thus, 
the  transmit  mode  uses  a  reflectarray  as  the  secondary  reflector  of  a 
two- reflector  system  and  the  receive  mode  uses  a  lens  with  a  simple 
reflecting  system.  The  basic  premise  on  which  this  system  is  based  is 
that  the  diode  phase  shifters,  which  comprise  the  phasing  control  of  the 
secondary  phasor,  can  be  made  to  operate  in  both  the  reflectarray  and 
lens  modes.  This  is  shown  schematically  in  Figure  2 

In  the  lens  mode  of  operation  (receive)  the  phasor  is  seen 
as  a  three-bit  balanced  hybrid  transmission  device  providing  45,  90, 
and  180-degree  increments  in  phase.  The  reflectarray  mode  (transmit) 
is  obtained  by  properly  biasing  the  180  degree  bit- -that  is,  one  diode 
is  forward  biased  and  the  other  is  reverse  biased.  Under  these  condi¬ 
tions,  the  180-degree  bit  becomes  a  short  circuit  and  reflects  all  the 
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FIGURE  2.  TRANSMISSION-REFLECTION  PHASOR 


energy  back  out  of  the  input  port.  This  provides  a  two-bit  phasor  having 
phase  increments  of  90  and  180  degrees  since  the  energy  passes  each 
bit  twice.  In  addition  to  providing  mode  switching,  the  180-degree  bit 
also  serves  as  a  duplexer  between  the  transmit  and  receive  modes. 

On  the  basis  of  this  simplified  description,  an  engineering 
study  is  being  performed  to  determine  the  design  parameters  for  a 
system  with  the  following  general  electrical  characteristics: 

TRANSMIT  MODE 


Frequency 

Half-power  beamwidth 
Zoom  factor 
Side- lobe  level 


3  Gc  +5  percent 
1x1  degree 
8 

-25  db  without  zoom 
-20  db  with  maximum  zoom 
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RECEIVE  MODE 

Frequency  3  Gc  +5  percent 

Half- power  beamwidth  1x1  degree 

Scan  angle  +8  degrees 

Side-lobe  level  -25  db  os  boresite 

-20  db  at  maximum  scan 
angle 

In  addition,  the  system  should  zoom  mid  scan  by  electronic 
means  and  be  feasible  for  implementation  at  UHF. 

The  effort  during  the  first  half  of  this  contract  has  been 
directed  toward  the  writing  of  a  generalized  computer  program  for  com¬ 
puting  the  far-field  radiation  patterns  in  three  dimensions  of  the  hybrid 
system.  This  program  will  serve  to  define  the  necessary  system 
geometry  needed  to  obtain  the  electrical  characteristics  stated  and  to 
determine  the  number  of  bits  needed  in  the  secondary  phasor  to  obtain 
the  zooming,  scanning,  and  side-lobe  level. 


SECTION  n 
SYSTEM  ANALYSIS 


Analysis  of  the  hybrid  system  ean  be  separated  into  two 
parts-  -determination  of  how  energy  propagates  within  the  system,  and 
calculation  of  the  far -field  radiation  patterns  using  the  results  of  the 
intrasystem  analysis.  This  division  is  possible  because  the  general 
characteristics  of  the  relationship  between  an  aperture  excitation  and 
the  resulting  far-field  pattern  are  fairly  well-known. 

An  experienced  observer  can  relate  the  two  functions  suf¬ 
ficiently  well  for  all  preliminary  design  and  the  more  refined  calcula¬ 
tions  can  be  used  for  a  final  detailed  analysis.  The  alternative  proce¬ 
dure  of  determining  the  scanning  and  zooming  patterns  for  each  change 
in  system  geometry  would  be  extremely  time  consuming. 

The  first  quarterly  report  of  this  program  contained  a  dis¬ 
cussion  of  the  intrasystem  analysis  and  the  computer  programs  needed 
to  perform  it.  This  report  '.ill  be  concerned  primarily  with  the  prob¬ 
lem  of  obtaining  the  far-fisld  radiation  characteristic  corresponding  to 
an  illumination  function  on  a  planar  aperture— the  output  from  the 
intrasystem  analysis. 

The  relationship  between  the  pattern  function  and  the  illum¬ 
ination  function  is  given  by  the  Fraunhofer  integral 


where 


x,y 
F(x,  y) 
u,v 


G<u,v)=  f  J 
A 


F(x,  y)exp(iux)exp(ivy)dxdy 


coordinates  of  planar  aperture 
illumination  function 

space  variables  2n  sin  9  sin  0/\  and  2irsin  9  cos  0/X 


6,0  =  spherical  coordinate  angles, 

Gfe,  v)  =  far-neid  pattern  function 

The  obliquity  factor  has  sot  been  included.  When  F{s,  y)  is  a  separable 
function  of  the  coordinates  the  pattern  function  becomes  the  product  of 
two  integrals  of  the  form, 

b 

r 

j  £(x)  esgi  {iux)  dx 
a 

in  inis  case,  the  problem  is  not  very  difficult  because  each  single  in¬ 
tegral  can  be  solved  rapidly  and  economically.  When  the  illuminaticn 
function  is  not  separable,  the  integration  is  difficult  and  time  consum¬ 
ing. 

One  approach  is  to  require  feat  the  antenna  systems  aper¬ 
ture  have  a  separable  illumination  function  under  all  conditions.  How¬ 
ever,  in  complicated  geometrical  configurations,  such  as  the  hybrid 
system,  the  energy  flow  is  quite  complex  and  such  a  constraint  would 
be  unrealistic,  leading  to  an  unfair  evaluation  of  fee  system's  capa¬ 
bilities. 

Most  attempts  to  obtain  a  general  method  of  numerically 
solving  fee  Fraunhofer  integral  have  approximated  fee  planar  aperture 
with  a  planar  array.  Solution  of  fee  array  factor  is  then  straight¬ 
forward  and  fee  precision  of  the  results  is  dependent  on  the  number 
of  array  elements  postulated.  Evaluation  of  the  total  array  factor  is 
quite  costly.  Performance  can  be  evaluated  in  several  planes  econ¬ 
omically,  but  intermediate  points  can  be  conveniently  solved  only  by 
changing  the  array  grid  structure  or  by  introducing  another  approxi¬ 
mation  with  an  undetermined  effect  on  precision  of  the  results. 


A  more  satisfactory  a sd  useful  solution  is  found  in  Hlsn’s 
method  (reference  3,  4)  of  quadrature  of  definite  integrals  of  the  form 


b 

|  f(x}  sin  udx  and 
j 

a 


cos  usds 


3 


sr 


Hie  method  involves  dividing  the  interval  (a,  b)  into  a  number  oi  equal 
parts  and  evaluating  sums  of  products  of  trigonometric  functions  and 
f(xc),  where  xQ  is  one  of  the  points  defined  in  the  interval  (a,  b).  The 
method  is  suggestive  of  tfcai  of  evaluating  array  factors  bur  the  com¬ 
puter  time  needed  is  markedly  less. 

Successive  application  ci  Rian's  quadrature  method  leads 
to  an  ajmrcsimaig  solution  of  the  Fraunhofer  integral.  A  more  com¬ 
plete  discussion  of  this  method  of  solving  the  double  integral  is  included 
in  the  Appendix.  The  precision  of  the  results  from  this  method  are 
primarily  determined  by  the  cube  of  the  spacing  between  grid  points 
in  the  aperture  and  by  the  fourth  derivative  of  the  illumination  function. 
Since  the  real  and  the  imaginary  part  of  the  illumination  func don  must 
be  evaluated  separately;  this  means  that  precision  degrades  approxi¬ 
mately  as  the  fourth  power  ef  the  scan  angle  and  as  the  third  power  of 
the  grid  spacing.  Best  results  are  obtained  when  the  aperture  is  sym¬ 
metrically  disposed  about  the  system  origin  (x  =  y  =  o). 

A  computer  program  was  prepared  (Appendix)  for  evaluat¬ 
ing  the  double  integral  by  Filon's  technique.  The  program  provided 
excellent  boresight  patterns  that  checked  exactly  with  known  results 
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for  uniform,  cosine,  and  cosine- squared  2pertare  illtmiinatioii  func¬ 
tions  and  for  combinations  of  these  functions,  for  example,  constant  x 
ees ins -squared.  The  accuracy  degraded  rapidly  with  scan  angle. 

The  results  of  a  study  of  the  spacing  required  to  obtain 
good  precision  for  a  diagonal  (worst  case)  scan  capability  of  13  degrees 
(0  =  45  degrees,  9  =  13  degrees)  are  shewn  in  Figures  3,  4,  and  5. 
Thirteen  degrees  was  considered  adequate  for  the  program.  An  il¬ 
lumination  function  of  cosine  x  cosine -squared  on  a  square  aperture 
was  used.  The  illustrations  show  the  principal  plane  patterns  for  the 
boresight  condition  and  for  the  diagonal  scan  condition  for  spacing  cf 
1,  1. 5,  and  2  wavelengths.  For  fee  least  spacing,  the  patterns  are 
almost  identical,  as  they  should  be  since  the  obliquity  factor  was  not 
considered.  As  the  spacing  increases,  fhe  scanned  pattern  deteri¬ 
orates  until  it  loses  almost  all  usefulness  at  the  maximum  spacing 
of  two  wavelengths.  The  computer  printout  for  one  wavelength  spac¬ 
ing  is  shown  in  Figure  6. 

A  pair  of  principal  plane  patterns  for  boresight  and  for  a 
26-degree  diagonal  scan  with  half  wavelength  grid  spacing  is  shown 
in  Figure  7.  The  agreement  is  excellent.  For  scan  angles  above 
30  degrees,  pattern  quality  degrades  rapidly  for  any  spacing. 

Considering  the  product  of  points  in  the  aperture  and  the 
far  field  to  be  Z,  the  present  cost  of  these  calculations  using  the 
IBM  1620  facility  at  AIL  is  about  Z/4  cents.  Improved  programming 
techniques  and  the  use  of  a  high-speed  computer  should  cut  this  figure 
by  more  than  an  order  of  magnitude,  probably  by  a  factor  of  50.  So 
the  present  cost  of  an  excellent  quality  13-degree  scanned  pattern  for 
a  60  x  60  wavelength  array  would  be  $9. 00/fer-field  point  with  improve¬ 
ment  to  less  than  $0. 25/far- field  point  obtainable.  The  cost  for  bore¬ 
sight  conditions  would  be  about  25  percent  of  the  above. 
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PRINCIPAL  PLANE  PATTERNS  FOR  TWO- WAVE¬ 
LENGTH  SPACING,  APERTURE  OF  20  x  20  WAVE- 
LENGTHS 
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SECTION  m 
CONCLUSIONS 


Analytic  techniques  for  investigating  the  performance  of 
the  hybrid  antenna  systems  have  been  studied  and  two  computer  pro¬ 
grams  have  been  prepared  for  calculating  the  essential  relationships 
between  the  antenna  configuration,  the  antenna  aperture  illumination, 
and  the  far-field  radiation  pattern.  Both  programs  have  been  tested 
against  known  results  and  have  shown  good  agreement  with  them. 
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SECTION  IV 

PROGRAM  FOR  NEXT  INTERVAL 

During  the  nest  quarter  of  this  study,  additional  hybrid 
antenna  configurations  will  be  analyzed  in  light  of  their  ability  to  pro 
•vide  the  required  pattern  characteristics. 

The  computer  programs  that  have  been  developed  will  be 
refined  and  put  in  proper  form  for  conjoined  use. 
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APPENDIX 

SOLUTION  OF  THE  FRAUNHOFER  INTEGRAL 
USING  FXLON’S  METHOD 

The  aperture  integral 

d  b 

G<u,  v)  =  f  f  F(x,  y)exp(iux)exp(ivy)dxdy 

can  be  numerically  solved  by  successive  applications  of  Filon's  quadra¬ 
ture  method  of  solving  oscillatory  integrals  of  the  form 


U  V 

f  2(x)  cos  uxdx  and  j  f(x)  sin  uxdx 
ar  a 

where  f(x)  is  continuous  in  the  interval  (a,  b). 

The  approximate  solutions  of  these  integrals  are  given  by 


if 

a 

where 


w  STS*-® 
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r  sin 


"V 
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cos  ub 
sin  ub 
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^r-l  =  Z  f(x2r-l] 
r=l  '  1 


cos  ux2r-l 
sin  ux2r_1 
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The  interval  (a,b)  has  been  divided  into  2m  segments,  each  of  length  h; 
the  segments  are  bounded  by  the  points  withs^  =  a 

and  x2rn  =  h. 


sin  8 
2Q2 


where  0  =  uh.  Expansions  of  these  functions  are  used  for  small  values 
of  0. 

For  the  double  integration,  the  interval  (a,b)  is  divided  into 
2n  segments  cf  length  h  and  the  interval  (c,  d)  is  divided  into  2m  segments 
of  length  l .  The  integral  must  be  solved  twice,  once  for  the  real  part 
of  F(x,  y)  and  once  for  the  imaginary  part.  Referring  to  either  part  as 
Q{x,  y)  we  have  the  double  integral 


which  leads  to  the  single  integral 


df  i 

h  J |ax |Q(a, y)  (-sin  ua  +  icos  ua)  +  Q(b,  y) (sin  ub  -  i  cos  ub) 
-jjp  J^Q(a,  y)  (cos  ua  +  i  sin  ua)  +  Q(b,  y)  (cos  ub  +  i  sin  ub)  j  + 


u 

Q  p2r*yj  [cos  uS2r  * 1  sin  ^r]  +  ~ 


yxL  ^(s2r-.P5j(dGSU3%r-l  +  isin^r-l)  /  *&&& 


This  expression  moot  be  integrated  term  by  term  with  respect  to  y.  An 
example  is  given: 

d  n 

%  f  y  Qp2r»y  cos  ux^_  exp(ivy)dy  = 
c  rd)  ' 


n  d 


SXcosux2rZ  f  «|x2r--y  exp(ivy)dy 


r=0  e 


n  f  r 

IPX  cos  uxgy  ffy  Q  c  (-sin  ve  +  i  cos  vc) 

r=0^  L 


Q  Xq  .  d  sin  vd  -  i  cos  vd 


Q|x2r»  c  (cos  vc  +  i  sin  vc)  +  Q  X2r,d  (cos  vd  +  isinvd)  + 
331 

$y  I  Q  |x2r’  *28  [ cos  +  1  sin  ’^s  ]  + 

o— A  \ 


m  f 

E  Q|x2r’y2s-l]  [cos  ^Zs-l  +  1  sin  ^s-l)  ) 


*yi.n 


Applying  the  same  techniques  to  the  oiherierms,  we  rert 
espress  the  double  integral  as 
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a 
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d  b 

J  QfeyJe^ir^eOTCivyJdxdy 


c  s 


where 


h+h*  +I11+I14  +  ]E16  +  —  +  124  + 
i{h  +  h  +— * I12  +  *13  +  X15  +—+  W 


u. 

X1  =  “^x  sinua  ^  Q(sj  y)  cos  vydy  s  -htc^  sin  ua(Vj) 
c 

d 

Ig  =  -haxsinua  J  Q(a,y)  sin  vydy  s  htc^  sin  ua(V2) 


d 

h  ~  8iavi>  J  Q(P>  7)  cos  vydy=  hlo^  sin  ub(Vg) 


u 

I4  =  hax  sin  ub  J  Q(b,  y)  sin  vydy  s  hio^  sin  ub(V4) 


d 

h  -  hf*x  L  cos  "^r  /  Qt*2r>*> 


cos  vydy  s 


W*  cos  uXgr(V5) 
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V  1 
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in  r 

IS  =  £  C03  nxgr  /  y>  sin  vydy  = 

r=Q  ^  ^ 


Ttl 


**x&'a"a*Fi 


i?_  f  in,p, 

I^-  =  -  ~2~  cos na  I  Q(a, y}  cos  vydy  s  — cos  iia(T/?) 


j 

c 


Ig  =  ~  -g*  cos  11a  J  Q(a,  y)  sin  yydy  =  — ^  cos  n&{Yg) 
c 


I8  =  - 


ljS_  ? 

-g—  cos  ub  /  Q(b,y)  cos  vydy  = 
c 


h£  r  nop 

ho  =  "  *ir'cos  ub  /  Qfe  y)  sin  vydy=  -  - -g  ■  cos  nb(VjQ) 


- TT^  COS 


nb(Vg) 


UR 
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X11  =  hyx  £  cosux2r-l  J  Q(x2r-l,y^  cos  s 

r=l  c 
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hi.yx  C03 
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■  f 

X12  =h^- £1C'0S  ^r-l  J 


m 

^xL  co?®2r-l(%) 


*13  =  fi£cts  cos  ■ca(V^) 


XX4=Mc££  003 


x15  =  cos  ^>(^3) 


X16  =  ~UaK  003  ^^4) 


hi  =  hwx  E0  ata  >“2^5) 


I18=htSxLsin'1Ji2r<V 


ht<3 

Ij!  9  =  -  — 2“^  Sin  ua(V7) 


hl%. 

X20  =  ~T'~sinua(V8) 


hi* 

*21  =  -  -~^siiiijb(VQ) 


UHX 

I22=”“2~sinub(vio) 
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ft 


A-6 


m 

*33  =  hvx  E  sin^r-lC’ll5 
r=i 

m 

%  *  hyx  E  Sin  “^r-l^ia5 

r=l 

For  each  far  field  point  these  24  terms  must  be  evaluated  and  properly 
summed  for  the  real  part  and  for  the  imaginary  part  of  the  illumination 
function-  These  results  are  combined  to  give  the  field  intensity  at  the 
point- 

A  printout  of  the  computer  program  used  to  perform  these 
operations  is  shown  in  Figure  8.  A  flow  diagram  of  the  program  is 
showD  in  Figure  9.  Improvements  in  the  programming  will  have  to  be 
implemented  before  the  program  will  operate  near  maximum  efficiency. 
When  this  is  done  the  integration  program  can  be  used  in  conjunction 
with  the  program  for  analyzing  the  energy  flow  within  the  system  in 
order  to  obtain  the  total  performance  characteristic  of  any  hybrid  sys¬ 
tem. 
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FIU3K  INTEGRATION  HH  TWO  DIMENSIONS  0010 

j*  AND  S  MUST  BE  COD 

DIHcHSIOH  X(12»12)»Y (12*12) tQ(12»12) »SH{30)  * SHI (301 ,CH{30 J  »CN1 (30) 
2*VKl30l *Fl£LO(3°i»COY(3P) 

R“A0  10fin»H,N»U0,UL.DU#V0*Vt»DV  0030 

OPAO  lfini.ysCAW.VSCAN 

READ  1002»A»r<*C»D  0050 

1000  FORMAT  (2I2,6?8.3>  2010 

1501  FORMAT C2FS* 31  2011 

1002  FORMAT (4F7«3J  2012 

EH=M-i 

£N=H-1 

H=(S-AJ/«  0080 

Et>(D-CJ/FN  0090 

Y2S=A 

DO  52  K=l»N 

COY{K)=*COSFt3.14159*Y2S/[B-Al> 

52  YZS=Y2$+RL 
Y2S*A 

W  31  K=!*N 
X2  S*C 

00  20  J=l.» 

CY=  3 . 1 4 1 54*C0SF ( USCAN*X  2S+VSCAN* Y2S ) 

S Y=-3. 14 1 59«S I NF ( USCAN»X2S+VSCAN*Y2S  J 
X  i  J  *  K ) =COY ( J  J  *CGY ( K 1 *C0 Y ( X 1 *CY 
Y  ( J  »K  >  =COY  ( J )  "•COY  (X 1  »COY  (X  1  *SY 
29  X2SaX2S+H 

31  v2Sxy2s+pL 
MlsN+1 
H1=M+1 

DO  301  X=1»N1 
X  (Ml  »*C  5  CQ*0 
301  Y<M1,K)»0.0 
DO  302  J=1 »M1 
X(J»NU«0,0 
30i  YJJ.N1 )a0.0 

U  (MK3  0100 

HF"(VL-V0)/DV+1.  0330 

V«VO 

0013  K“1*NF 
VK(X)«V 
13  VsV+OV 

PRINT  1010»(VWXJoK*l»NFl 
1910  FORMAT  (9X,2HV*,13F8.3I 
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FIGURE  8.  PRINTOUT  OF  FELON  2  PROGRAM  US  TING 
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KF=lUL-UOJ/fHJ+l. 

DO  100  I=V,HF 

14  TH=(U*HJ 
TTH=2i«TH 
THSQ»TH*TH 
THCU=THSO*TH 
TKR  =S£J  RTF!  THSO  1 
IF£THS>-.25115»16*16 

15  ACX=(.04444444-(. 006349206-. OG04232804»THS01*THSOJ*THCU 
BEX*. 6666667+1. 1333333-1 .03809524-.003527337*TH5Q>*THSQ)*THSQ 
GAX=1 .333333-1 .1333333-1 »004761905-«0OQ08818342*TH5QJ#THSQ>*THSO 
60  To  17 

16  CTTH=COSr  1  TTHl 
STTH=SINF!TTH* 

CTH*C0SF1TH) 

STH=SINF(TH> 

AJ_X“  ( 1  ,+STTH/TTH-l  1  ,-CTTH) /THSO) /TH 
3EX=  ( 3  ,+C TTH-2.  +STTH/TH  >  /THSO 
6AX=-4.* I CTH-STH/THJ /THSO 

17  UA=U*A 
UB=U*B 
V=VO 

DO  90  L=1 *NF 

18  PH=iV*FU 
TPH=2.«PH 
PHSO=PH*PH 
PHCU=PHSO»PH 
PHR=SQR  TF ( PHSO ) 

20  !F(PHR-.25 )21 »22 *22 

21  ALY=< ,04444444-i .006349206-. 0004232804*PHSQ)*PHSQ»*PHCU 
BEY=. 6666667+1. 1333333-{.03809524-.003527337*PHSQ>*PH.SQ>*PHSQ 
GAY=1 • 333333-1 .133333 3- ( . 004761 90 5-. 000086 18 342*PHSO ) *PHSO)*PHSO 
GO  TO  23 

22  CTPH=COSF ( TPH 1 
STPH=SINF1TPH) 

CPH»COSF(PH) 

SPHbSINF(PH) 

ALY* ( l.+STPH/TPH- 1 1 ,-CTPH  > /PHSO ) /PH 
BEY*13.+CVPH-2.#STPH/PH)/PHSn 
GAY«-4.# ( CPH-SPH/PH > /PHSO 

23  VC«V*C 
V0-V*0 

24  S1-5INF1UAJ 
S2=SINF1UB) 
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S3»St«HVC) 

SA-SfUFCVDl 

C1*C0SRUA1 

.C2“COSc(U3j 

C3*C05FIVCS 

C4eC0SFJ¥Dl 

S7»53 

S8*S4 

C7“C3 

C8«=G» 

25  Y2S=A 
Y2S1=A+FL 

26  DO  3«  K*1,N,2 
K1=K+1 

SN(K)=S!Nr (V*  Y2SJ 
SfJlOCl  )=SINF(V*Y2S1  J 
CN  ( K ) =COSF ( V*Y2S  J 
Cm  (KlMCOSr  CV*Y2S11 
Y2S=Y2S+2.#FL 
30  Y2S1=Y2S+FL 
DO  70  J= 1 *  Ml  - 
DO  70  K=1»N1 
70  0(J.KI=X<J,K) 

LN=0 

SUMX=0. 

SUV.Y = 0  « 

32  Rl=-.5*{0(l,l)*C7+onvN)»C8l 
R11=0» 

R2=-.5*(0(M,l>»C7+0(w,N)*C8) 

P21=o. 

ra=-, 5* ( 0 ( 1 , 1 } *C7+0( 1 ,N ) #C8 ) 
R41=n. 

R5e-.5»(0(K»l)*C7+0(MtN)*C8) 

R51-9. 

Tl=-.5*(OU,l)*S7+0(l  ,N)*SB) 
T11=0, 

T2«-.5*(0tM*l)*$7+0(M,N)*S8) 

T21-0. 

T4=-.5*(0(1.1>*57+0(1»N1»S8> 

T41=i. 

T5*-,5»{0(M.l ’s$7+0(M»N)*S8) 
T31«0 

34  DO  40  K»ltN»2 
K1=K+1 
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Pl»9t+OU»KJ»CW(KJ 

Rn*pii4oci,ia)»c»ncKT  j 

1010 

m 
«’ ■& 

R2»R?+0(H.K)«CN(K) 

1020 

R21=R21+0«MtKH«CNl(KU 

1030 

R4*R4+0  *  I *X 1 *CN  f K ) 

1040 

R41*P4l+Q(l,>a  )*CN1(K1 i 

1050 

14 

U 

R5*B540tM,K)*CN(KI 

1050 

R$laR5l40(M,Kl)*CNl (K1 1 

1070 

T1«T140(1,<)»SM(SC) 

1030 

Vi 

Tll=TH4(}n*K'*  )*SN1(K12 

1090 

«5  * 

*-< 

T2=T2+0(M.K)*;,N(K) 

1100 

T21=T2l40iM,KX)*SNl(tCll 

1110 

h 

T4=T4+0{ 1*K)*SN(K) 

T41=T41+0(1»K1)*SN1(K]  > 

1130 

pi 

T6=T6+0(M,K,>#SN(K1 

1140 

40 

T51=T51+0(M,K1 )*SN1(<1 ) 

1150 

\ip% 

42 

HL=H*FL*ALX 

1160 

H 

V1=ALY* (-0(1*1) *S340 ( 1 ,N ) *S4 ) 4BEY*R 1+GAY«R1 1 

1170 

K  a 

FI1=-V1*HL*S1 

1180 

f V 

V2=ALY* (Q( 1  *  1 7*C7-Q( 1 *N)»C8 )+BEY*Tl+GAY*Tll 

1190 

4?. 

EI2=-V2*HL*S1 

1200 

f-J 

V3=ALY*(-Q(M,1)*S7+Q(M,N)*S8)+BEY*R2+GAy*R21 

1210 

i*j 

FT3=V3*HL*S2 

1220 

r  .* 

V4=ALY*(Q(M,1 !*C7-0(M*N)*C8)+BFY*T2+GAY*T21 

FI4=V4*HL»52 

1240 

0 

N 

46 

HLn=H»FU»BEX 

1250 

V7=AlY» (-0(1*1 ) *S3+Q{ 1 ,N )*S4 ) +BEY«R4+GAy*R41 

1270 

EI?=-V7*HL3*»  5*0 

1280 

<  y 
'/  it- 
! 

V8=ALY* ( 0( l , 1' *C3-0( 1 *N  >  *C4 ) +BEY*T4+GAy*T41 

1290 

E18*-V8*HLD*.:»C1 

1300 

v  5 

VO~AUY*(-0(K,ll*S3+0fM,N)*S4)+BEY*R5+GAY*R5l 

1310 

i 

F-I0=-V«)*HLB*,A*C2 

1320 

jmJ 

Vin-ALY*(0(M,1)*C3-0(M,N)*C4)+BFY*T5+GAY*T51 

F! 10=-V10*HLB*«5*C2 

1340 

El  1 3<=Vl*HL*Cl 

1350 

El  14=-V2*HL*C1 

1360 

\&C. 

El  15s-V3*HL*C2 

1370 

FI  16=V4*HL*C2 

1380 

i 

FI14«-V7«HLB*.5*S1 

1390 

F120bV8*HLP*.5*S1 

1400 

P!21«-VO*huB».5*S2 

1410 

E122*  VI 0*HLP*»5*S? 

1420 

■^0- 

48 

FI5-0 

1430 
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ft  <y 


A-12 


r  c*  ^  V 


E111*0 

PH2«0 

.FU?«0 

EM8-0 

£123*0 

£124*0 

X2R*C 

X2R1*'>K 

DO  SO  J»1#P*2 

JI*J+1 

R3*-i 5»  C  0 1 J . 1 > *C7+0{ J  ffN ) *C8 ) 

R31*0 

R6=-<,5*f0{Jl»ll*C7+0{Jl.N)*C8) 

R61*0 

T3=-«5*<0( J»I J*S7+OC JjN)*S81 
T31*0 

T6*-.5otO(Jl.- :  J*S7+r(Jl*N)«S8) 

T61=0 

DO  PO  K*1»N.? 

Xl-X+1 

R3=R3*Q<J.KJ*CNtKl 
R31=R31+OC J.K1 >*CN1 1  tC  > 

R6=R6+Q(J1.K»*CN<K> 

R61=R61+0(Jl»KlJ*CNmit 

T3=T3+0(J»K**SN<K) 

T31=T31+OtJ*Kl )*SN1 t XI S 

T6=T6+0(JlfK)*SN<»U 

SO  T61=T61+Q( J1 »K1 > #SPil  UO  i 

V5=ALY*  t-0( J. 1 ) *S3+0( J.N ) «S4 ) +BEY*R3+GAY*R31 
EJ5=EIP+JHLB*V5*COSFtU*X2R) J 
V6=ALV*<Q(J»1>*C3-Q{ J,N>*C8 >+BEY*T3+GAY*T31 
FI6*Et  6+Hi.B#V6*COSF  ( U*X2R  > 

HLG=H*EL*GAX 

Vll=ALY*  ( -Q<  J1  «  1 )«S7+0( J1  »N  I  *S8  J  +BEY#R6+GAY*R61 
E!U*EI11+HL0*VU«C0SFCU»X2R1) 
VX2oALY»(r>tJX»lJ#C7-OlJl,N)*C8)  +SFY#T6+C-AY*T61 
6112*FI  12+HLG'*V12*C0$F{U*X2R1  > 
EIX7«FtX7+KLB<-V5*SINF<U*X2R) 
E!18*FT18-HU3»V6*SINF  «U*X2R5 
E!23«E!23+HLG*V11  *SINF(U»X2R1 > 
EI24«FT24-HIG*V12»STNF(U*X2R1> 

X2R«X2R+2.»H 
60  X2R1 *X2R4H 
rF(LN)61. 61.62 
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61  SUMX=FI  1+FI3+E1 5+EI7+FIP+EI 11+FI 14+FI 16-frE! 18+FI20+E1 24+EI22  1870 

SUMYsEI 2+EI4+EI6+F18+E! 10+EI 12+EI 13+EI 15+EI17+EI 19+FI21+EI23  1880 

LN=1  1890 

DO  71  J=1 »M1 

DO  71  K=1 sNl 
71  OIJ,K)=Y( J,K) 

GO  TO  32  1900 

62  SUMX=SU1'X-E 1 2-E 1 4-E 1 6-E 1 8-F II 0-E I 12-EI 13-EI 1 5-E 1 1 7-F 1 1 9-EI 21-El 23 

SUMY=SUMY+EI 1+FI3+EI3+FI7+EI9+FI11+FI14+EI16+EI 18+EI20+FI 22+E 1 24  1920 

30OO  FORMAT  (2F8.3) 

F ! ELD { U  =4. *SQRTF  (  SUMX*SUMX+SUiV.Y*SUMY  ) 

QO  V=V+DV  2020 

PRINT  l005,U9(FIELD!l.)»L=l*NF) 

1003  FORMAT  (/3H  U= ,F3 .2 . 1 3F8 . 3 ) 

100  U=U+DU  2030 

64  CONTINUE 
CALL  EXIT 
FND 
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